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Abstract 

We show that dimensionful renormalization scheme parameters such as 
the renormalization or factorization scale can be completely eliminated from 
perturbative QCD predictions provided that all the ultraviolet logarithms 
involving the physical energy scale Q are completely resummed. 
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The problem of the renormalization scale (scheme) dependence of fixed- 
order perturbative QCD predictions continues to frustrate attempts to make 
reliable determinations of the underlying dimensional transmutation param- 
eter of the theory , Aq C d (usually or a s (M z ) are the fitted quantities). 
Whilst a number of proposals for controlling or avoiding this difficulty have 
been advanced [I], 0, 0, f| no consensus has been reached, with the result 
that in experimental fits attempts are made to estimate an ad hoc "renor- 
malization scale" uncertainty |J . This "theoretical error" can be larger than 
the actual experimental errors, and in our view can potentially mislead as 
to both the central value of and the likely importance of uncalculated 
higher-order corrections @, || . 

It is undeniable that with standard fixed-order renormalization group 
(RG-) improvement there is a scheme dependence problem. Any proposed 
"solution" must, therefore , amount to special pleading for a particular choice 
of scheme, motivated by albeit reasonable considerations imported from out- 
side perturbative field theory. Examples include the BLM approach of Brod- 
sky and collaborators f3[ in which the piece of the next-to-leading coefficient 
proportional to the first beta-function coefficient is absorbed into the cou- 
pling, motivated by various QED examples; or the Principle of Minimal Sensi- 
tivity (PMS) criterion of Stevenson [ED] , where the scheme is to be chosen such 
that the perturbative approximation is as insensitive as possible to changes 
in scheme. In this paper we wish to emphasise that the renormalization scale 
dependence of fixed-order QCD perturbation theory is due to the incomplete 
nature of the standard RG- improvement carried out. The idea will be that 
the dependence of a dimensionless QCD observable 7Z{Q) on the dimension- 
ful parameter Q (e.g. the cm. energy in e + e~ annihilation) is obviously 
completely independent of how the theory is renormalized. The perturbative 
coefficients and the coupling a s (fj,) in contrast manifestly depend on the di- 
mensionful renormalization scale fi , via the presence of logarithms ln(/i/A) 
, where A is universal and depends on the subtraction procedure used to 
absorb infinities (usually MS). The Q-dependence, however , is built by ul- 
traviolet (UV) logarithms ln(Q/A^) with A-ji completely independent of the 
renormalization procedure. The asymptotic Q-dependence of 7Z{Q) is deter- 
mined by the dimensionful parameter An which is a physical property of the 
observable, independent of the renormalization scheme @. The key obser- 
vation is that one should keep \i independent of Q. If this is done standard 
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fixed-order RG-improved predictions are self-evidently inadequate since they 
do not satisfy asymptotic freedom , 1Z(Q)— >0 as Q— >oo. This property only 
results if all the RG-predictable ln(Q/An) terms are resummed to all-orders 
building a l/ln(Q/An) behaviour. This summation may be accomplished 
with any choice of /i , but in doing so all /i-dependence cancels between the 
ln(/i/A) and a s (/j,) terms. Thus the complete resummation of all the RG- 
predictable UV logarithms gives /^-independent predictions. One has traded 
unphysical /i-dependence for the correct physical Q-dependence. 

In the standard RG- improvement one truncates the resummation of the 
In^Q/An) terms and uses a Q-dependent scale /i = xQ, this "incomplete" 
improvement yields an x-dependent result. Clearly the UV logarithms are 
physical and leaving out an infinite subset of them gives the resulting fixed- 
order approximation an unphysical scale dependence. Of course there is still 
residual scheme dependence since fixed-order predictions will depend on the 
other dimensionless parameters specifying the scheme, which can be taken to 
be the non-universal beta-function coefficients JIJ , but dimensionful scheme 
dependence parameters can always be eliminated by complete resummation 
of all Q-dependent logarithms. 

There are close links with the Effective Charge approach of Grunberg [0 
which focusses on building the Q-dependence of TZ(Q) . Incomplete improve- 
ment with fi chosen to be the effective charge (or fastest apparent conver- 
gence (FAC)) scale is equivalent to the "complete" RG- improvement outlined 
above. Similar remarks apply to examples such as moments of leptoproduc- 
tion structure functions where there are two (or more) dimensionful scales, 
for instance the renormalization scale \x and in addition a factorization scale 
M. In this case one has ln(M/A), ln(/t/A) logarithms as well as "physical" 
UV logarithms involving Q, which are independent of the renormalization 
and factorization conventions. Again all /t and M dependence is eliminated 
provided that all the physical UV logarithms are resummed [0] . 

We begin by briefly reviewing the problem of parametrizing RS-dependence, 
and define the concept of RG-predictable terms. Consider the dimensionless 
QCD observable TZ(Q) , dependent on the single energy scale Q (we assume 
massless quarks). Without loss of generality, by raising to a power and scal- 
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ing, we can arrange that 7Z(Q) has a perturbation series of the form, 



Tl(Q) =a + na 2 + r 2 a 3 + . . . + r n a n+1 + . . . , (1) 

where a = a s (p)/n is the RG-improved coupling. The /i-dependence of a is 
governed by the beta-function equation, 

dct 

-ba 2 (l + ca + c 2 a 2 + . . . + c n a n + . . .) . (2) 



<9ln/i 

Here h and c are the first two universal terms of the QCD beta-function 

33 - 2N f 

b = — —> 

153 - 19N f 

r = — 

126 



(3) 
(4) 



with Nf the number of active quark flavours. As demonstrated by Steven- 
son the renormalization scheme may be completely labelled by the vari- 
ables r=61n(/i/A) and the non-universal beta-function coefficients C3, . . .. 
a(r, c 2 , c 3 , . . .) is obtained as the solution of the transcendental equation 

\ + cln {iT7a) =T - [ dx (-^ + ^(iT^)) ' (5) 

where B (x) = x 2 {l + cx + c 2 x 2 + c 3 x 3 + ...). This is obtained by integrating 
Eq. (2) with a suitable choice of boundary condition related to the definition 
of A. 



For our purposes it will be useful to label the RS using n, the next-to- 
leading order (NLO) perturbative coefficient, rather than r. This is possible 
because 

r- ri = Po (Q)=b\n(Q/A n ), (6) 

where po is an RS-invariant, hence r can be traded for r x . is a dimension- 
ful scale dependent on the particular observable. It is related to the universal 
dimensional transmutation parameter A^r by 

A^ r/i A® , (7) 
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where r=rj us (/i = Q). The righthand side of Eq.(7) is independent of 
the subtraction scheme employed, and as advertised has a physical sig- 
nificance, being directly related to the asymptotic Q-dependence of H(Q) 
[§, H, 0. Equations (5), (6) can then be used to define a(r 1 ,c 2 ,c 3 , . . .). 

We next turn to the RS-dependence of the perturbative coefficients ?v 
This must be such as to cancel the RS-dependence of 'a' when the series 
is summed to all-orders. The self-consistency of perturbation theory |J de- 
mands that the result of a N n LO calculation (terms up to and including 
r n a n+1 ) in two different schemes should differ by 0(a n+2 ). This implies 
the following dependences of the r« on the scheme parameters- rzirijC?), 
r 3( r i, C2, C3),. . ., r n (ri, c 2 , C3, . . ., c n ). The self-consistency requirement can be 
used to derive expressions for the partial derivatives of the r n with respect 
to the scheme parameters. For instance for r 2 one has 

dr 2 dr 2 dr 2 , , 

— = 2rx + c, — = -1, — = 0, ... . 8 

OT\ OC 2 OC 3 

On integrating these expressions one finds 

r 2 (n, c 2 ) = r x 2 + cr 1 + X 2 - c 2 

5 1 

r 3 (ri, ca, ca) = r x 3 + -cn 2 + (3X 2 - 2c 2 )n + X 3 - -c 3 

; ; . (9) 

In general the structure is 

r n {ri, c 2 , . . ., Cn) = r n (n, c 2 , . . ., c n _i) + X n - c n / (n - 1) . (10) 

Here f n is an n th order polynomial in r\ which is determined given a complete 
N n_1 LO calculation. X n is a Q-independent and RS-invariant constant of 
integration and can only be determined given a complete N n LO calculation. 
f n is the "RG-predictable" part of r n , and X n is "RG-unpredictable" . Thus, 
given a NNLO calculation in the MS scheme with \x = Q one can determine 
the RS-invariant 

X 2 = rp(^ = Q) - (rp(/x = Q)f - c r™QjL = Q) + cf^ , (11) 

where 

77139 - 15099iV f + 325iV 2 

Co = / • (12) 

2 17286 V ; 
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By a complete N n LO calculation we mean that c 2 , C3, . . ., c n have been com- 
puted as well as ri, r 2 , . . ., r n . We note in passing that c^ s has now been 
calculated ||. 

Using Eqs. (9) we can now exhibit the explicit RS-dependence of the terms 
ofEq.(l), 

5 1 
TZ(Q) = a+r 1 a 2 +(r 2 1 +cr 1 +X 2 -c 2 )a 3 +(rl+-cr 2 1 + (3X 2 -2c 2 )r 1 +X 3 --c 3 )a 4 +. . . 

(13) 

where a=a(ri, c 2 , c 3 , . . .). We adopt the principle that at any given order of 
Feynman diagram calculation all known (RG-predictable) terms should be 
resummed to all-orders. Given a NLO calculation r x is known but X 2 , X 3 , . . . 
are unknown. Thus the complete subset of known terms in Eq.(13) at NLO 
is 

5 1 

a =a + r x a 2 + (rf + cr x - c 2 )a 3 + {r\ + -cr\ - 2c 2 r 1 - -c 3 )a 4 + (14) 

The sum of these terms, ao , can be simply determined using the follow- 
ing two-step argument. The infinite subset of terms in Eq.(14) has an RS- 
independent sum, since the X 2 , X 3 , . . ., -dependent terms cannot cancel their 
RS-dependence, and we know that the full sum of Eq.(13) is RS-invariant. 
Each term is a multinomial in ri, c 2 , C3, . . .. Using the RS-independence we 
can set r\ = 0,c 2 = 0,c% = 0, . . ., in which case all terms but the first in 
Eq.(14) vanish and we obtain ao = a{r\ — 0, c 2 — 0, C3 = 0, . . ., c n = 0, . . .). 
So at NLO this corresponds to working in an '"t Hooft scheme" with c 2 = 
c 3 = . . . = PT0|] , and with r\ = 0. From Eq.(6) r x = corresponds to 
r = b\n(Q/ An) or to an MS scale [i = e~ r ^ b Q. This is the so-called "fastest 
apparent convergence" (FAC) or effective charge (EC) scale [I], |2[. From 
Eq.(5) we find that ao satisfies 

l. + c]n(-^-)=b]n(Pi . (15) 
a Vl + cao/ \An/ 

We note to avoid confusion that the definition of A on which Eq.(5) is based 
[|]] differs from that usually used for A^r. In terms of the standard definition 
we have, 

, /?r\ ~ c l h 

A^ = e^ -) A— s . (16) 
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If a NNLO calculation has been completed , then X 2 can be determined 
(as in Eq.(ll)), and a further infinite subset of terms are known and can be 
resummed to all-orders, 

X 2 a 3 = X 2 a 3 + 3X 2 na 4 + . . . . (17) 

The RS-independence of the sum and the multinomial structure of the co- 
efficients again leads to a resummed result involving a . We finally arrive 
at 

K(Q) =a + X 2 a 3 + X 3 a 4 + ...+ X n a n+1 + ..., (18) 

which is simply the perturbation series in the RS with r\ = c 2 = C3 = . . . = 
c n — . . . — 0. As is obvious from Eqs.(9), X n = r n {r\ = 0, c 2 = 0, . . ., c n = 0). 



Unfortunately the result obtained by resumming all RG-predictable terms 
depends on our choice of r±, c 2 , C3, . . ., c n , . . . as the parameters used to la- 
bel the scheme. Whilst this choice is natural and straightforward it is evi- 
dently not unique. We could equally consider a translated set of parameters-: 
f\ = r\ —r±, c 2 = c 2 — c 2 , . . ., c n = c n — c n , where the barred quantities 
are constants. The partial derivatives in Eq.(8) with respect to these new 
parameters are unchanged so that 

p = 2f I+ 2r 1 + c, p = -l, p = 0, (19) 
or 1 oc 2 0C3 

which on integration yields 

r 2 (r 1 ,c 2 ) = r\ + 2r 1 r 1 + cf 1 + X 2 - c 2 

; ; (20) 

with general structure 

r n {r x ,c 2 , . . .,c n ) = f n (f 1 ,c 2 , . . .,c n _i) + X n - c n /(n - 1) . (21) 

The X n are again constants of integration which are unknown unless a com- 
plete N n LO calculation has been performed. If one applies the same rationale 
as before, where all RG-predictable terms are to be resummed, one finds , 
analogous to Eq.(18), 

K(Q) =a + r 1 a 2 + X 2 a 3 + ... + X n a n+1 + ... , (22) 
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with a=a(r\ = Fi, c 2 = c 2 , ■ ■ ■, c n = c n , . . .), which is just the perturbation 
series in the scheme with r\ = r±, c 2 = c 2 , . . ., c n = c„, . . ., or equivalently 
with fi = 0, c 2 = 0, . . ., c n = 0, . . .. 



Thus by itself the principle of resumming all the RG-predictable terms at 
any given order of Feynman diagram calculation does not abolish the scheme 
dependence problem since the subset of RG-predictable terms depends on 
how the scheme dependence is parametrized, the parametrization ambiguity 
being precisely equivalent to the scheme dependence ambiguity. However, as 
we shall now argue, the parameter r\ has a special status, being connected 
with the dimensionful renormalization scale fi and the physical energy scale 
Q on which TZ(Q) depends. Thus from Eq.(6) we see that we may write 

rM= (bin J- -bln-f] , (23) 



with /i taken to be the MS scale as is customary. To simplify the discussion 
let us temporarily set c = and work in a scheme with c 2 = C3 = . . . = c n = 
. . . = 0. Then the coupling a(/z) is given by 

«w=v^wy, (24) 

and the NLO RG- improvement in Eq.(14) becomes a geometrical progression 
in n, 

1l(Q)na(n) + n(ii)a 2 (n) + r\{^)a\^) + ...+ r> n+1 (^) + . . . . (25) 

Substituting Eq.(23) for ri(/i), summing the geometrical progression and 
using Eq.(24) for a(/i) yields 



H(Q)^a(n)f 



1 — f & In -J^- Mn ) a(a) 



l/61n(Q/A w ) . (26) 



We see explicitly the cancellation of the unphysical ln(fi/Aj^g) terms in r 1 ( / u) 
with those in a(/i) to generate the correct physical Q-dependence of 7Z(Q), 

1Z(Q)nl/b\n(Q/A n ) + 0(l/61n(Q/A^)) 3 . (27) 

In contrast with standard NLO RG- improvement one has 

n(Q)^a(fi)+ ri (fi)a 2 (fi) , (28) 
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where \i = xQ is taken proportional to Q with x a dimensionless constant, 
and x = 1 the so-called "physical scale" is often favoured. The resulting 
Q-dependence is 



KiQ^l/biHxQ/Ajfs) +0(l/b\n(xQ/A W s)f , 



(29) 



which is, of course, x-dependent. All x-dependence cancels and the physical 
Q-dependence of Eq.(27) results if the geometric progression of Eq.(25) is 
not truncated but is resummed to all-orders. Truncating the resummation 
can result in considerable error in the extraction of Ajjs from comparisons 
of NLO perturbative results with experiment. For many e + e~ jet observ- 
ables one has r = rf IS (fi = Q)~10 Q and the truncation of a geometric 
progression with a common ratio ra~l/2 (taking a = = Mz)~0.05) 

can lead to a sizeable overestimate of the true unless fortuitously the, 
as yet unknown , NNLO and higher invariants X2-,X^ 1 . . . compensate. 

The complete resummation of Eq.(25) is forced on one if the renormal- 
ization scale /i is kept independent of the physical energy scale Q (i.e. /i 
held constant) since then the ln(Q/A^) UV logarithms are the only source 
of Q-dependence. Standard fixed-order NLO RG-improvement is then man- 
ifestly unsatisfactory since it does not satisfy asymptotic freedom, 1l(Q)—>0 
as Q^oo. Instead 



tends to — oo as Q— >oo with /i fixed. Only if the complete resummation in 
Eq.(25) is carried out does one build the correct 1 /ln(Q / A-ji) behaviour. 

Notice that use of the parametrization with r± = r\ — Ti and T\ = 
b\n(xQ / Ajjs) — bln(Q/Aji) is equivalent to standard NLO RG-improvement 
with fi = xQ and yields the x-dependent result in Eq.(29) when all the NLO 
RG-predictable f i terms are resummed. The resummation is then manifestly 
incomplete with respect to UV logarithms since the ria 2 term is split into 
fia 2 + Tia 2 with the r^a 2 term not summed. Furthermore the constants of 
integration X 2 ,X 3 , . . ., which are unknown and hence omitted at NLO now 



R-NLoiQ) = a{fj) + bin T 




(30) 
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contain \n(xQ / 'Ajjg) and ln(<5/A^) terms. For instance 

X 2 = X 2 +(b\n^--b\n^-) . (31) 



If these terms are included in the resummation the x-dependence disappears 
order-by-order reproducing the complete resummation of Eq. (25) and giving 
the correct physical Q-dependence in Eq.(27). So the parameter t\ has a spe- 
cial status and should not be translated if one wishes to completely resum 
all the UV logarithms and reconstruct the correct physical Q-dependence. 

With realistic non-zero c the NLO complete RG-improvement (CORGI) 
of Eq.(14) sums to ao satisfying Eq.(15) . This may be written in closed form 

as 



K(Q) 



•1 



c[l + W(z(Q))\ 



<Q) = --(^-\ bl \ (32) 



e \A n 

where W(z) is the Lambert iy-function defined implicitly by W(z)exp(W(z)) ■ 
z [|J. This is of course equivalent to standard "incomplete" NLO RG- 
improvement with jx = e~ r / b Q, the FAC (or EC) scale , as noted ear- 
lier. Whilst CORGI yields //-independent results there is still a depen- 
dence on the other dimensionless scheme parameters c 2 , C3, . . .. That is one 
should parametrize the RS-dependence using r 1 (i.e. r 1 = 0), but there 
is no preference for any particular parametrization £2,03, .... In the effec- 
tive charge approach of Grunberg one chooses c 2 ,C3,...,c n so that the 
X 2 , X 3 , . . ., X n are zero at N n LO, corresponding to r\ — 0, r 2 — 0, . . ., r n = 0. 



The use of c 2 , C3, . . ., c n as parameters corresponds to an 't Hooft scheme |L0 
c 2 = 0, c 3 = 0, . . ., c n = 0, and is a priori equally reasonable 0. At NLO one 
conventionally makes this choice in any case. 

We finally note that inverting the relation of Eq.(32) leads to the property 
of asymptotic scaling [O] 



lim QF{1l{Q)) = A n , (33) 

Q^oo 
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where T{x) is the universal QCD scaling function 

F(x)=e- 1/bx (l + l/cx) c/b . (34) 

Exactly this property is used in lattice gauge theory calculations of lattice 
coupling to assess how close to the continuum limit of infinite inverse lattice 
spacing one is. By using the exact relation between A-ji and the universal 
parameter Ajjg in Eq.(16) one arrives at the universal asymptotic scaling 
relation 

lim QFCJl(Q))e- r/b (2c/by /b = A^ . (35) 

Q— >oo 

This relation can be used analogously to the lattice scaling relation to assess 
how close to asymptotia in Q one is at current energies, for QCD observ- 
ables calculated to NLO. One simply inserts the data for TZ(Q) , and the 
corresponding NLO corrections r=r^ s (fi = Q) , into Eq.(35). The asymp- 
totic prediction is that all the data should lie on a single horizontal straight 
line corresponding to Ar^g. Deviations from this unambiguously indicate the 
presence of sub-asymptotic effects , and enable the estimation of relative dif- 
ferences between the uncalculated NNLO invariants X 2 , and possible power 
corrections, for different observables. Observation of the scaling property 
at some approximate level provides a well-founded starting point for further 
analyses attempting to resum large logarithmic corrections for jet observ- 
ables or predict power corrections. In particular the sub-asymptotic effects 
are contained in a non-universal factor Q{TZ(Q)) which approaches unity as 
Q^oo, so that 

Qf(nQ)MnQ))e- r/b (2c/b) c/b = A^ , (36) 

with 

g{K{Q)) = 1 - (X 2 /b)K(Q) + ... + (Ko/n^K-^e- 1 /™ + .... (37) 

The K term represents a possible power correction, here taken to be 1/Q. 
Direct fits of X 2 and K to the Q-dependence of the data can be performed 
along the lines discussed in Ref.[6]. 

This is to be contrasted with standard experimental analyses which at- 
tempt to assess the size of uncalculated higher-order corrections by variation 
of the chosen renormalization scale /i = xQ. As we have attempted to indi- 
cate here such an approach can be misleading , and leads to no information 
on the size of the uncalculated higher-order invariants X 2 , X%, .... 



11 



Acknowledgements 

We would like to thank John Campbell and Nigel Glover for a productive col- 
laboration and numerous discussions which were instrumental in developing 
the ideas reported here. 



12 



References 

[1] P.M. Stevenson, Phys. Rev. D23 (1981) 2916. 

[2] G. Grunberg, Phys. Lett. B95 (1980) 70; Phys. Rev. D29 (1984) 2315. 

[3] S.J. Brodsky , G.P.Lepage and PB. Mackenzie, Phys. Rev. D28 (1983) 
228. 

[4] D.T. Barclay, C.J. Maxwell and M.T. Reader, Phys. Rev. D49 (1994) 
3480. 

[5] S. Bethke, Nucl. Phys. B (Proc. Suppl.) 64 (1998) 54. 

[6] J.M. Campbell, E.W.N. Glover and C.J. Maxwell, Phys. Rev. Lett. 81 
(1998) 1568. 

[7] C.J. Maxwell and A. Mirjalili, in preparation. 

[8] Einan Gardi and Georges Grunberg, JHEP 9807 007 (3776328) (1998). 

[9] S.A. Larin, T. Van Ritbergen and J. A.M. Vermaseren, Phys. Lett. B400 
(1997) 379. 

[10] G.'t Hooft in : "The Whys of Subnuclear Physics", Erice (1977), ed. A. 
Zichichi (Plenum, New York, 1977). 

[11] C.J. Maxwell, Phys. Lett. B409 (1997) 450. 



13 



